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Production of Drell-Yan pairs and open and hidden heavy flavor on nuclear targets 
is examined within perturbative QCD. The effects of modifications of nucleon structure 
functions inside the nuclear medium are considered. Besides, nuclear dependence of 
charmonium and bottonium absorption is studied in the framework of the Glauber- 
Gribov model. The low energy limit of this approach recovers the probabilistic formula 
usually employed for charmonium and bottonium suppression in nuclear collisions. 



^^Talk given at the 5th International Workshop on Hard Probes in Nuclear Collisions (Lisbon, 
Portugal, September 7th- 18th 1997). 
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1 Introduction 



Drell-Yan (DY) and open (QQ) and hidden (if?) heavy flavor production in hadron col- 
lisions are usually studied in the framework of perturbative QCD [[l], Q]. The extension 
to nuclear projectiles and targets is straightforward, provided factorization between 
partonic densities and parton-parton cross sections is assumed: 

f 1 f 1 

OAB = <J AB ^ hH = J2 dx * dx b fa/A(Xa^ 2 )fb/B(x b ,/J 2 ) a ab ^ hh (§, m h , // 2 ) . (1) 

a b JXaO Jx b0 

In this equation the summation runs over all partons in the projectile and target, f c /o 
is the density of parton c in hadron or nucleus C, a ah ~* hh is the parton-parton cross 
section and factorization and renormalization scales are taken to be equal; = rriQ 
(M/+/-/2) for heavy flavor (Drell-Yan) production, x a0 = 4ml/ s, xjo — 4m|/(sx a ), 
s = sx a Xb and /i 2 ~ m 2 . The leading order (LO) relation xf = x a — x b will be used. 

In this context, some differences between A, B being hadrons or nuclei are the 
following: 

• The influence of the nuclear medium on nucleon structure functions or partonic 
densities (see for example 0), expressed as the ratio: 

R^) = AUid(x ^I- (2) 

If we parametrize cross sections on nuclei as <Jab(xf) — &pp (AB) a ( XF \ this effect 
makes a(xp) ^ 1 for all processes, even if isospin effects (i.e. the difference 
between neutron and proton parton densities) are corrected. 

• The scattering of the produced heavy system (absorption by nuclear matter, 
rescattering with hadronic co-movers and/or deconfinement, see 0), which affects 
if? production, making '(xp) < a DY (xp). 

• The elastic scattering of initial partons, which is the accepted explanation of the 
Pr-broadening ((pt)ab greater than (pt) pp for DY and \P production), being the 
difference proportional to A 1 ^ ||, 
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• The energy loss of fast partons inside the nuclear medium (the so-called jet 
quenching), proportional to A 2 ^ 3 and different for cold and hot nuclear matter 
II 0]- It may affect the yield of charm and bottom at high energies |§. 

In this contribution we will examine the two first aspects. It is organized as follows: 
In Section ||] nuclear structure functions will be briefly discussed, and the parametriza- 
tions used in our calculations and results for Drell-Yan, open heavy flavor and charmo- 
nium production on nuclear targets will be presented. In Section |3] nuclear absorption 
of states with hidden heavy flavor will be studied. Finally, in Section |] some conclusions 
will be presented. 



2 Hard processes on nuclear targets 

Since the experiments of the European Muon Collaboration ||, the modification of 
nucleon structure functions inside the nuclear medium has been studied by several 
experiments ||. While the dependence of R a /A(x, Ai 2 ) on /i 2 is very small, four regions in 
x can be described: R a /A > 1 for 0.8 < x < A (Fermi motion and cumulative regions); 
R a/A < 1 for 0.3 < x < 0.8 (the original EMC effect); R a/A > 1 for 0.1 < x < 0.3 
(antishadowing region); and R a /A < 1 for x < 0.1 (shadowing region). 

We will use the parametrization of Ref. |10]]Qfor Au, which was designed to describe 
the four regions in x. This approach follows the conventional techniques for global fits 
to produce nucleon parton distributions: the ratio R a / A is parametrized at some low 
virtuality (4 GeV 2 ) and then evolved to higher /i 2 using evolution equations modified for 
the nuclear case [jn|. All parameters are fixed from a comparison to nuclear structure 



function ratios over deuterium and DY data, using baryon number and momentum 

sum rules and a SU(3) symmetric sea. This parametrization has a lower limit in 

x = 10~ 3 (= ^/s ~ 100 GeV for charm at xf = 0, not enough for predictions for RHIC 
2 An update of this parametrization can be found in jflj; among other modifications, modern sets 
of nucleon structure functions have been used. Other proposal in the same spirit can be found in |B3 . 
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and LHC or for high xp = x a — x^), so at smaller x we have taken two alternatives: 
either frozen or a linear-log extrapolation in the form x". The parametrization for two 
different virtualities can be seen in Fig. 1 [p~4| . 

Figure 1: Parton densities in Au for fi 2 = 5 and 25 GeV 2 : valence quarks (dashed 
curve), sea quarks (solid curves) and gluons (dashed-dotted curves); dotted curves are 
the results of linear-log extrapolations with (3 = 0.096 and 0.040 at jj 2 = 5 and 25 
GeV 2 respectively. 
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2.1 Open heavy flavor and Drell-Yan 

Parametrizing ct p a{xf) = v PP A a ( XF \ in Figs. 2-5 results |l4j] for a(xp) in DY, c and 
6 production in pAu collisions are presented, using GRV HO |15j and MRS A |16| 
nucleon structure functions. The a ~* (s, m^, /i 2 ) are taken at next-to-leading order 
in the MS renormalization scheme and the following parameters are used: 

i) For open heavy flavor |I7|| : m c =1.5 GeV and /i 2 = 4 GeV 2 for charm, and m& = 5 
GeV and /i 2 = m 2 for bottom. 

ii) For Drell-Yan |TJ: /j, 2 = M 2 +l _. 
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As stated above, Xp = x a —Xb (LO relation), x a xt,s > 4m| and the main contribution 
to the integrals comes from x a XbS = 4mf, x a /b = [\J {lQm\/ s) + xp ±Xf)/2] this means 
x a /b — 3 • 10~ 4 for charm at yfs = 10 TeV and xf = 0. 

Figure 2: Energy dependence of a for charm and beauty production in pAu collisions 
for GRV HO (solid and dashed curves) and MRS A (dotted and dashed-dotted curves) 
structure functions and using extrapolated (dashed and dashed-dotted curves) and 
frozen at x = (solid and dotted curves) ratios of parton distributions. 

cx 

1 

0.98 
0.96 
0.94 
0.92 

2 3 4 

10 10 10 

Vs NN (GeV) 




It can be seen in these Figures that the influence of the chosen Set of nuclear parton 
densities is almost negligible except for the highest energies and xp. Besides, some 
difference appears between frozen and extrapolated ratios; as expected, this difference 
is larger at higher energies or xp and for charm or low dilepton masses. Moreover, in 
the xp distributions all regions in nuclear structure functions can be seen from negative 
to positive xf (corresponding to decreasing x). Related results can be found in Refs. 

m 
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Figure 3: xp dependence of a for charm and beauty production in pAu collisions at 
= 39 GeV (upper figure) and 1800 GeV (lower figure) for GRV HO and MRS 
A structure functions and using extrapolated and frozen at x = 10~ 3 ratios of parton 
distributions (with the same conventions as in Fig. 2). 
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2.2 Hidden heavy flavor 

We will concentrate on charmonium production on nuclear targets. Now the for- 
mation of the final resonance has to be considered. Usually two models are used to 
study charmonium production: 

i) The Color Evaporation Model (CEM) |BlJ considers that all color dynamics is con- 
tained in the kinematical restriction to Eq. (fj): Am? c < s = x a XbS < Am 2 D . The 
projection on different charmonium states is taken into account by numerical coeffi- 
cients, ff AB ^ = a AB ^ cc (4ml < s < 4m|,), which are universal in this approach. 
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Figure 4: Mass (upper figure) and energy (lower figure) dependence of a for Drell-Yan 
pair production in pAu collisions for GRV HO and MRS A structure functions and 
using extrapolated and frozen at x — ratios of parton distributions (with the 

same conventions as in Fig. 2). In the lower figure M 2 is in GeV 2 . 
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ii) The Factorization Approach (FA), which contains both the Color Singlet Model 
(CSM) and Color Octet Model (COM), is based on non-relativistic QCD (NRQCD) 
H |22|. In this model the parton-parton cross section for production of a charmonium 
state ^ is: 

a"(J,m c , / i 2 )=E^ W (Oj ] ), (3) 
M 

with C^Li the short distance coefficients for the hard subprocess ab — > cc[n] ([n] is 
the color configuration of the cc pair), computable as series in a s , and (0*j) the long 
distance matrix elements taking into account the hadronization cc[n] which can 
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Figure 5: xp dependence of a values for heavy lepton pair production in pAu collisions 
at ^snn = 39 GeV (upper curves in each figure) and 1800 GeV (lower curves in each 
figure) and different masses for GRV HO and MRS A structure functions and using 
extrapolated and frozen at x — 10~ 3 ratios of parton distributions (with the same 
conventions as in Fig. 2). Note that at ^s^n — 39 GeV all curves coincide. 
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be classified in powers of the relative velocity v of Q and Q and are obtained from a 
fit to data [23, 24, 25 1. We will use the fit to fixed target data obtained in Ref. |2^ . 



While in CEM different resonances have the same energy and Xp behavior, in FA 
we consider different contributions and take into account the most important decays 
from charmonium states into J/'0 : 

a t J % = af^ + B(^^J/^)a^+ £ B{ XcJ — J/4>) <x XcJ , (4) 

J=0,l,2 

being a d jj^ the direct J/ip production, i.e. not coming from decays, and B(H — > J/ip) 
the branching ratios for particle H to decay into a J/ip. Experimentally it is found 
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that direct J/ip contribution is about 60 % of the total J/ip cross section, if)' decay 
gives less than 10 % and all Xcj contribute with more than 30 %. Decays into if>' are 
not important, and then if)' production is dominated by direct production. 

One important point is the contribution of the different color states to the produc- 
tion of these particles, i.e. the color content of the pre-resonant state. In fact, the FA 
gives that direct J/if> production is almost completely produced in color octet state, 
if}' is also predominantly (about 90 % or more) produced in color octet, and the main 
contribution of the Xcj states to J j if> comes from color singlet matrix elements. Then, 
a separate study of the effect of nuclear structure functions on different particles and 
color states is possible in this approach. As we will see, the color octet and color singlet 
contributions to the production of charmonium have different suppression. 

In Figs. 6 and 7 we present results for energy and xp dependence of charmonium 
production p6| . Computations have been done taking m c = 1.5 GeV and \i = 2m c , 
and using CTEQ3L [E7[] and GRV HO [TjJ nucleon parton distributions in FA and 
CEM respectively. Different color contributions and charmonium states are taken into 
account defining 

af(x F ) = (5) 

for i = CSM, COM, CSM+COM=FA, CEM, and ^ = total J/ip, direct J/ip, if)', 
J2B(xcj — ► J/i') a xcj- ft can be observed that both FA and CEM give very similar 
results except for x's at very high energies or xp. The behavior of the different color 
contributions (and hence of the S wave and P wave states) is also different at high 
energies or xf- 



3 Nuclear absorption 

Nuclear absorption of the pre-resonant cc state in its path through nuclear matter 
is usually taken into account in pA collisions, at fixed impact parameter b, by two 
formulae (which neglect nuclear effects on structure functions, assumed very small at 
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Figure 6: Center of mass energy dependence of a for pAu collisions. Different lines are: 
FA total contribution (solid), singlet contribution (dotted), octet contribution (dashed) 
and CEM (dashed-dotted). 




fixed target energies in the central rapidity region): 



i) Taking into account the path L(b) across the nucleus | 2q| , Fig. 8 



(Tp A (b) oc exp[-p a abs L(b)}, p ~ 0.17 fm" 



(6) 



ii) The Glauber probabilistic formula f29fl : 

(jp N A / d 2 b j dz p(z, b) exp 

J J— oo 



pA 



-OabsA 



dz'p(z', b) 



&abs 



d% 



1-e 



-o a hsAT A (b) 



(7) 



this formula can be easily understood (Fig. 8): the pre-resonant cc state is created 
at some point z in the nucleus (in an amount proportional to the density of nuclear 
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Figure 7: xf dependence of a for y/s = 39 GeV pAu collisions with the same conven- 
tions for lines as in Fig. 6. 




matter p(z, b) at that point) and is absorbed in its path through the nucleus from z to 
+oo with some absorption cross section^ For open charm {a abs = 0), ct p a = <7 p n A. 

Eqs. @ and fl7|), with their longitudinal ordering, are only valid in the low en- 
ergy limit. A formula valid at all energies has been derived |HJ in a fully relativistic 
Glauber-Gribov approach, using finite energy Abramovsky-Gribov-Kancheli (AGK) 

cutting rules [55] . Considering n nucleon-nucleon interactions taken in some arbitrary 
3 The meaning of this absorption cross section is not clear. In some proposals (e.g. pof ) it has been 
related to the color structure of the pre-resonant state: absorption is much stronger in octet than 
in singlet configuration. Observation of little absorption for x' s ( see Subsection 2.2) and variation of 
absorption with px of the produced resonance in some defined way would support this point of view 
and the validity of the NRQCD approach. 
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Figure 8: Representation of the absorption mechanism at low energies in the rest frame 
of the nucleus; z is the creation point of the pre-resonant cc pair. 




longitudinal ordering zi < z 2 < ■ ■ ■ < z n , these rules imply (besides the usual AGK 
prescription [Q) to change T^(6), T^fe) = f_™dz p(z,b), by 



/+oo r+oo r+oc 

dzi j dz 2 -- ■ / dz n exp [zA(zi - Zj)\ \\ p(z i: 6), (8) 
-oo J z\ Jz n -i 

for j the first interaction either cut or to the right of the cut (see Fig. 9) and 

(9) 



m N M 2 , , 
A = , M = 2mh or m^,. 



s.i 



This is nothing but the t m i n effect ||34j| , i.e. the nuclear form factor suppresses interac- 
tions which require a minimum momentum transfer (e.g. those in which heavy flavor 
is produced). 

Both an external contribution (corresponding to the heavy system being produced 
in the interaction) and internal^ contributions (corresponding to the heavy system 
already present in the projectile or target, which turn out to be very small) are included. 
Modifications of the nucleon structure functions due to the nuclear medium naturally 
appear in this framework, which is valid at not very high x. The results for the low 
and high energy asymptotic limits are: 



J Not to be confused with intrinsic charm pa], which can be important at high x. 
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Figure 9: Position of the cut to the left of the j-th interaction (dotted line) or cut in 
the j-th interaction (crossed solid line). 



j-i 




i) For low energies, A — ► oo ==>- only j = 1 contributes, the shadowing of structure 
functions disappears and we recover the probabilistic formula (|7|). To our knowledge 
it is the first time that this equation has been derived in a fully relativistic approach. 

ii) For s — > oo, A = and T^\b) = T^(b), the shadowing of the structure functions 
factorizes and we get for the external part: 



pA 



a 



a 



d 2 b e 



-uAT A {b)l2 



I _ e -«AT A (b)/2 



a* N A / d 2 b T A {b) e 



-aAT A (b)/2 



(if no shadowing). 



(10) 
(11) 



a (<t) is the light (heavy) particle- nucleon cross section and <j a bs = (1 — e)cf (e can 
be interpreted as the probability for the heavy particle to survive in one interaction, 
e = 1 (~ 0) for open charm (charmonium)). Terms with a correspond to the mod- 
ification (shadowing) of nucleon structure functions inside nuclei; in Eq. (|10]) it has 
been described by an eikonal model using a multi-pomeron factorized vertex, although 
other models could be used (as a sum of fan diagrams, i.e. the Schwimmer model, see 



mm- 

The consequences of this approach, neglecting the modifications of nucleon structure 
functions inside nuclei, are the following: at high energies open charm and charmonium 
are equally absorbed (no a a b s appears in Eqs. ( |T0"D and ([□])), see also [0); the low 
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(0) and high (|TT|) energy formulae differ for charmonium up to 20 % at the highest 
energies (being the difference of order (&A/R\) 3 , a ~ o" a 6s); and the exact [[H]] an d 
probabilistic results differ ~ 1 4- 2 % at = 20 GeV. 

These results can be interpreted as follows: At low energies only the first interaction 
is effective for producing heavy flavor (so production is proportional to A), subsequent 
interactions can only absorb it. At high energies all interactions are simultaneous, so 
the absorption mechanism of subsequent interactions is no longer effective; instead the 
full multiple interaction formalism, which suppresses equally both hidden and open 
heavy flavor production, has to be considered. Usual factorization (i.e. separation 
between partonic densities and parton-parton cross sections, Eq. ([!])) is broken and an 
additional suppression factor is always present (except for open heavy flavor production 
at low energies). While at finite energies it is not possible to separate in this additional 
factor the modification of nucleon parton densities inside nuclei from the scattering of 
the heavy partons, this separation is recovered at low energies (Eq. (Uf), where there is 
no modification of nucleon parton densities) and also, for fixed impact parameter b, at 
asymptotically high energies (Eq. (|IUD). 

Results for pPb collisions ]31] , in the form of the variation with energy of A e ff = A a , 



are presented in Figs. 10 and 11. We use a oc (xs) (normalized to 7 mb at pi 



ab 



200 GeV/c [§]]), e = 1 (0.001) for open (hidden) charm and a standard Woods-Saxon 
nuclear density. Just to give some estimation, the effect of nuclear structure functions 



has been taken into account considering it factorized and computed following |LJ, |26 
see Section ||. It can be observed that the exact result provides a smooth transition 
between the low energy and the high energy regimes (the latter already reached at 
yfs = 200 GeV) and that the effect of the nuclear modification of structure functions 
varies, at xf = 0, from antishadowing at low energies to shadowing at high energies. 

Turning back to the xp dependence of charmonium nuclear absorption, it is clear 
(see for example |ft7| ) that the modification of the nucleon parton distributions inside 
nuclei cannot reproduce the data: for different energies, there is no scaling of absorption 
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Figure 10: A e ff versus y/s for xf = for charmonium and open charm production in 
pPb collisions: exact result |H| (solid line), probabilistic formula (Eq. (7), dotted line), 
asymptotic formula (Eq. (11), dashed line) and exact result with nuclear modifications 
of structure functions as explained in the text (dashed-dotted line). 




in Xb (corresponding to the target nucleus), while there is approximate scaling in xf = 
x a — Xb- Nevertheless nuclear structure functions have to be taken into account. In 
Fig. 12 calculations [^(| of a versus Xp for pW collisions at pi a b = 800 GeV/c are 
presented, with the same parameters as in Subsection 2.2; normalization is obtained 
using Eq. (0) with 

Cabs — 6 mb (taken constant with xp). Clearly the modification 
of the nucleon parton distributions inside nuclei accounts for part of the effect, except 



at the highest xp. More data, both at 800 GeV/c, as those from [ 381 , and at lower 
energies, have to be examined (following e.g. the proposals of |31], |34], |39|) in order to 
explain the behavior of the absorption with xf- 

As a comment, let us stress that a is a misleading variable. For example, in terms 
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Figure 11: The same as Fig. 10 but for xp = 0.5. 
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of a the experimental results of [|37]] and seem to indicate some anomalous behavior 
at xp — in collisions at 800 GeV/c: absorption increases going from xp ~ 0.15 to 
xp ~ 0.65 and also from xf ~ 0.15 to xp ~ 0. If one expresses these results in terms of 
the ratio W/C, which can be done both for the E772 [[37]] and E789 |38j data, this ratio 
turns out to be monotonically decreasing with xp from xp ~ to 0.65 (0.780 ± 0.074, 
0.776 ± 0.059, 0.746 ± 0.046, 0.741 ± 0.020, 0.729 ± 0.034, 0.649 ± 0.051, 0.604 ± 0.084, 
0.571 ± 0.155 for x F = -0.023, 0.032, 0.16, 0.26, 0.36, 0.46, 0.55, 0.65 respectively, 
with only statistical errors taken into account). So, it follows the trend that can be 
seen in Fig. 12 and no strange behavior can be deduced. 



4 Conclusions 



Effects on heavy flavor and Drell-Yan [14] and charmonium |26|] production of the 
modification of nucleon structure functions inside nuclei have been examined. They 
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are relatively small at low xp, low yfs but become very important at high xp and for 
RHIC and LHC energies. Thus, more experimental and theoretical effort is needed 
to reduce the uncertainties in the extrapolation to small x. For charmonium nuclear 
absorption is, at low energies and xf ~ 0, of greater importance than nuclear effects 
on parton densities; for higher energies or larger xp both effects have to be taken into 
account. 

Figure 12: Comparison of xp dependence of nuclear corrections by modifications of 
parton densities inside nuclei with experimental data of 800 GeV/c protons incident 



on a tungsten target from [37]. Theoretical calculations have been normalized to one 
of the less shadowed experimental points (see text). Lines follow the same convention 
as in Fig. 6. 




Formulae for nuclear production of open and hidden heavy flavor have been pre- 



sented [31], using a relativistic Glauber-Gribov formalism in which the standard prob- 
abilistic formula for charmonium absorption has been derived as a low energy limit of 
an exact expression valid at all energies. The numerical accuracy of the probabilistic 
formula at available energies has been checked. A striking prediction of this approach 



31 , |34}| is that, at high energies, open charm will asymptotically be as suppressed as 



charmonium. 
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New data |4(J on J/ ip suppression in PbPb collisions at SPS energies have produced 
great excitement as a possible signal of new physics (i.e. Quark-Gluon Plasma). In 
view of all the uncertainties commented in this contribution, detailed tests of all our 
conventional ideas about open and hidden heavy flavor production off nuclei are needed 
before any quantitative statement on the existence of new physics in heavy ion collisions 
can be made. 
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